Controlled hole doping of a Mott insulator of ultracold fermionic atoms 
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Considering a system of ultracold atoms in an optical lattice, we propose a simple and robust im- 
plementation of a quantum simulator for the homogeneous t-J model with a well-controlled fraction 
of holes x. The proposed experiment can provide valuable insight into the physics of cuprate super- 
conductors. A similar scheme applied to bosons, moreover, allows one to investigate experimentally 
the subtle role of inhomogeneity when a system passes from one quantum phase to another. 
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Gases of ultracold atoms in optical lattice poten- 
tials provide extremely clean and controllable condi- 
tions for studying strongly correlated many-body physics 
Since for deep lattices these systems are described 
quantitatively by Hubbard- type Hamiltonians @, they 
have great potential to serve as quantum simulators [3[ 
for paradigmatic models of condensed matter physics. 
This prospect is fed by results for bosonic systems. 
In a seminal experiment, the quantum phase transi- 
tion from a superfluid of bosons to a strongly corre- 
lated Mott insulator — predicted for the bosonic Hubbard 
model [i[ — has been observed Moreover, quantitative 
agreement between experiment and ab initio quantum 
Monte Carlo simulations clearly confirm the validity of 
the Bose-Hubbard description @. 

While the elementary Hubbard model for bosons is 
rather well understood, this is not the case for repul- 
sively interacting fermions: In a Mott insulator, with 
interaction localizing one particle of "spin" s = f or J, 
at each site, Fermi statistics gives rise to an antiferro- 
magnetic superexchange coupling J between neighbor- 
ing spins. Intriguing physics is expected when such a 
quantum antifcrromagnet is frustrated, either by a non- 
bipartite lattice geometry Q or by doping it with holes 
(or fermions) displacing spins when moving around [Toj . 
The latter scenario for a square lattice is conjectured to 
give rise to d x 2_ y 2--w&ve pair superfluidity and to explain 
basic properties of high-temperature cuprate supercon- 
ductors [a4lfj|. However, conclusive theoretical evidence 
of whether the plain Hubbard physics supports a super- 
conducting state is still lacking. As pointed out already 
in Ref . [ll[ , here a cold atom realization of the fermionic 
Hubbard model (or its descendant, the t-J model) could 
provide critical insight. 

The recent observation of a Mott insulator of repul- 
sively interacting fermionic atoms in a deep optical lat- 
tice [l2j is an important first step toward a clean cold 
atom realization of strongly correlated fermionic Hub- 
bard physics. However, further steps in that direction 
require solutions to two problems: (i) The temperatures 
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that can be achieved presently are still larger than (or 
at most comparable to) the energy scale of the superex- 
change spin-spin coupling (l3| . Promising novel cooling 
procedures have been proposed to tackle this problem 
fl3 | (cf. also [Hj]). (ii) While it is relatively easy to cre- 
ate an incompressible Mott region with one atom per 
site in the center of a parabolic trap, it is very hard to 
dope such a trapped Mott insulator in a controlled way. 
When the particle number in the center of the trap is 
lowered, e.g. by slowly widening the trap, the central 
Mott-insulator phase will not be doped with holes ho- 
mogeneously. Rather, it will melt from the edge. Also 
simply switching off the trapping potential (i.e. compen- 
sating it with a blue-detuned laser) is difficult, since then 
particles may leak out and it will be hard to control the 
lattice filling. In this Rapid Communication we propose 
a robust solution to problem (ii) . Our method allows one 
to both accurately control the fraction x of hole doping 
and effectively compensate the trapping potential. Both 
are crucial in order to learn about high-temperature su- 
perconductivity with ultracold atoms. Applied to bosons, 
our scheme, moreover, allows one to experimentally in- 
vestigate the non-trivial influence of spatial inhomogene- 
ity when a system passes from one quantum phase to 
another [l6|. 

Our basic idea is to create holes in a fermionic Mott 
insulator by adding auxiliary bosonic particles to the sys- 
tem that interact repulsively both with the fermions and 
with each other. Each boson repels a fermion from one 
site. Experimentally, one has to create a mixed Mott 
phase with one particle, either fermion or boson, per site. 
This spatially rather well-defined "simulator region" re- 
sembles the system to be simulated. The degree of hole 
doping is well controlled by the number ratio between 
bosons and fermions. The price to be payed is that, un- 
like the tunneling of a real hole, the motion of the bosons 
happens via boson-fermion swaps, b eing slow second- 
order superexchange processes (cf. [Tj, [l8| for the case of 
"spinlcss" Bose- Fermi mixtures). But this type of hole 
kinetics brings also a great advantage: Consider the trap 
being very similar for fermions and bosons (think of two 
Ytterbium isotopes in a far off- resonant dipole potential) . 
Then such a boson-fermion swap on two neighboring sites 
will change the potential energy only marginally, even if 
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FIG. 1: (color online) (a) Sketch of the lattice and trapping 
potential felt by both fermions and bosons. Assuming strong 
on-site repulsion and symbolizing fermions by light red ar- 
rows and bosons by blue bullets, an occupation snap shot is 
sketched. There is a central mixed Mott region, the simu- 
lator region, with one boson or fermion per site, (b) In the 
simulator region, fermions move by fermion-boson swaps not 
changing potential energy. Here a description solely in terms 
of a new type of fermions (dark red arrows) is possible, with 
bosons implicitly taken care of as holes. These fermions move 
with an effective tunneling amplitude and do not feel the trap; 
the fraction x of holes equals that of bosons. 



the trapping potential does change between both sites. 
Hence, within the simulator region the physics is hardly 
influenced by the trap, cf. Fig. [TJ Below we will show 
that this region is accurately described by the homoge- 
neous t-J model with controlled hole fraction x. The t-J 
model [lH (cf. [20[ for other cold atom applications) de- 
scribes the doped fermionic Mott-insulator, with fermion 
hopping t and superexchange spin coupling J (« t/3 in 
cuprates [2l[ ) . Second-order superexchange is crucial for 
the physics of the doped fermionic Mott-insulator. We 
want to emphasize that our proposal, in which both spin 
coupling and hopping originate from superexchange, does 
not involve any process above second order. Our scheme 
is, in a sense, contrary to the slave-boson approach to the 
t-J model (where auxiliary bosons describing holes are 
introduced as a purely theoretical concept): We propose 
a physical system of real bosons and spin 1/2 fermions 
that behaves as a system of fermions alone. 

Before presenting the details of our proposal, we want 
to mention that recently another, interesting solution to 
problem (ii) of controlled hole doping [not to the issue of 
temperature, problem (i)] has been proposed [22]: The 
authors suggest realizing an attractive Hubbard model, 
which for a bipartite lattice can be mapped to a repulsive 
one (see, e.g., [23]). In particular, an imbalance of t 
and 4- fermions for attractive interaction, being relatively 
easy to control experimentally, corresponds to hole or 
particle doping in the repulsive model. Unfortunately, at 
the same time, the trap felt by the attractive fermions 
transforms to an inhomogeneous magnetic field, favoring 
a spatial separation of repulsive f and J, fermions. 

Let us now describe in detail our approach to controlled 
hole doping. We consider a mixture of Nr = Nf + N± 
fermionic and Nb bosonic atoms (N = Nf + Nb) in an op- 



tical lattice, with the fermions having equally populated 
"spin" states s =j, l. The system is described quantita- 
tively by the Hubbard model Q 

Bff* = -E [ E t iMu+ h - c -)+ t U h lh+ h - c -) 

jjbb 

+ Y [ uff ^it^i,l + U fb n. lf n. lb + —n t b(h ib - 1)] 

i 

+ J2 [Vifn + sVihif]. (i) 

i 

Here f} s and b\ denote creation operators for fermions 
and bosons at site i. We also define number opera- 
tors h is = fi s f is , nib = b\b v hif = h it + h^, and 
hi = fiif+hib. The first line of Eq. ([1} comprises fermion 
and boson tunneling between neighboring sites, with pos- 
itive matrix elements t{j and tjy (depending exponentially 
on lattice depth 2]). The second line takes care of the 
repulsive on-site interaction between the different types 
of particles, with Hubbard energies U ff , U fb , and U bb 
(weakly depending on lattice depth and being propor- 
tional to the corresponding s-wave scattering lengths (2|). 
Finally, the third line includes co-centric trapping poten- 
tials V/ = Vi+ 5Vi and V b = Vi for fermions and bosons 
[V./ = V b = in the center] , with 6V, to be tuned small. 

We are interested in the parameter regime giving rise 
to an extended mixed Mott region, with one particle 
(boson or fermion) per site, in the center of the trap. 
We assume a cubic lattice of spacing d with site loca- 
tions dri, as well as a parabolic confinement Vi — \arf. 
We also define both U min = mm(U ff , U fb , U bb ) and 
t max = max({ty}, {t\j}). Temperatures are well lower 
than i max . Now, i max < U min and p, < U min (while 
t max <C U m in — A 4 )) with chemical potential fx, guaran- 
tees strong suppression of double occupancy. On the 
other hand, the trap Vi prevents vacancies from enter- 
ing the central region: the occupied region of radius 
p ss (3./V/47T) 1 / 3 features only a thin shell (a few d wide) of 
reduced particle number, provided i max < ap = AT4d ge , 
where AV e d g e is the energy difference between sites at 
radius p and p + 1. With the above requirements ful- 
filled, the chemical potential for bosons and fermions is 
basically given by the potential energy needed to place 
a particle at the edge of the occupied region, p w \ a P 2 ■ 
One has AI4dgc/M = 2/0. The radii p might range from 
10 to 30. All in all, i max < AV^dge -C p < U m - m sum- 
marizes the parameter regime assumed here. In practice, 
imax/kmin can be adjusted via the lattice depths, while 
/V^min can be tuned to a value of 1/2, say, by varying 
the trap depth a. 

In the bulk of the mixed Mott region, at each site i, 
strong repulsion in combination with the trapping con- 
finement Vi gives rise to the constraint 

hi = h it + h. ti + h ib = 1, (2) 

with the overall boson fraction x ~ Nb/N. The system 
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can be described by an effective Hamiltonian H e ff act- 
ing in the subspace Si defined by Eq. ©. Treating the 
first line of Hf/b as perturbation Hi , we can expand H c s 
according to degenerate perturbation theory. 1 Up to sec- 
ond order, using fin, = 1 — hif in Si, one finds 

h cS = f{-^E M/W5/i. + h - c -) + E want 

(ij) s i 

+ Y: [M&iSi ^) + U™n lf n 3f ] }P, (3) 

(ij) 

where P projects on Si and where we have introduced the 
usual spin operators Si = hY^s's Ha ,(T s'shs with Pauli 
matrices oy s . Moreover, 
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with % = [(^) 2 - I# = 4%£[1 + 8\% 

rfb = (*{ 7 ) 2 + (^) 2 n c/fci „/6 = («f ; ) 2 -(^) 2 r/b ^ 
J ij — z [7/6 L 1 ~r" °jj J > 0J ij — [/7 b °ij Vt-Vj ' 

and adj(i) containing all sites adjacent to i In Hamilto- 
nian ([3]), tij describes boson- fermion swaps, i.e. effective 
fermion tunneling, Wi is the effective potential felt by 
fermions, Jij stands for the usual fermionic antiferromag- 
netic superexchange coupling, and Uf™ captures boson- 
mediated nearest-neighbor interaction between fermions 
that can be attractive, repulsive or zero. 

Hamiltonian (|3]) is equivalent to a t- J-type model, de- 
scribing a purely fermionic system when strong repulsion 
suppresses double occupancy. This can be seen by intro- 
ducing composite-fermion creation operators 

4 s =flk- (8) 

We define h is = c\ s c is and Si = \ J2 S > S ^\s lCT s'sC is . Using 

Pflf is P = PflA^ + VMsP = MsAsP, y^lds both 
His = hi s and Si = Si in Si . Bosons transform to empty 
sites (holes) and Si to the subspace defined by hi < 1 
with hi = hrf + hi±. Thus, we can rewrite if e ff as t-J 



1 With unperturbed states a, a' £ Si, b S\ and en- 
ergies E a , E a i, Eft, the leading orders of H e ff read [24)1 : 
(a'\H^\a) = (a'\H \a), (a'\H^\a) = (a'\Hi\a), and 
(a'\H^\a) = T,bW\Hl\b)(b\Hi\a)\[(E a - E b )~ 1 + (E a , - 
Eb)- 1 ]- 



Hamiltonian in terms of cJ s -fermions alone: 

HeS = P{-J2Y, *«(4c j - a + h-C.) + J2 W ^i 

(ij) s i 

+ ]T Jij [(SiSj + U^hihj] }P. (9) 

(ij) 

Transformation © (being inverse to a the slave-boson 
transformation) is illustrated in Fig. [JJ 

The t-J model © with W l = const., UTf = 0, U 3 = t, 
Jij = J « t/3, and boson fraction x, is the most simple 
candidate to explain high-temperature cuprate supercon- 
ductivity (§l-Hol. Hf| . The fermion-boson mixture consid- 
ered here, allows one to realize these parameters quite 
accurately: Starting from a cubic lattice, the tunneling 
amplitudes t(j and t\ - can be suppressed in one direction 
by ramping up the lattice in that direction, leading to a 
stack of uncoupled square lattices layers. From now on, 
we will only consider the intra-layer physics. Creating 
the optical lattice by using a rather broad laser beam, 
within the occupied region of the trap (created by fur- 
ther beams) the lattice will be practically homogeneous. 
Thus, t{j ~ tf and t\ 3 - ~ t b . Furthermore, potential differ- 
ences | Vi — V 3 ; | between neighboring sites are much smaller 
than U ff , U fb and U bb , giving 5£ < 10~ 3 (v = ff, fb, bb) 
for the parameters estimated in a previous paragraph. 
With that the model parameters ((4j)- ((6J are to good ap- 
proximation homogeneous within the simulator region, 
tij ~ t, J^ ~ J, and UVf ~ U nn . By the very same 
arguments, also the last three terms contributing to Wi 
[Eq. (Jll)] have a negligible spatial dependence compared 
to the effective hopping parameter t, being the relevant 
energy scale here. Thus, if also the difference between 
boson and fermion potential SVi is smaller than t, one 
achieves a practically flat effective potential Wi ~ W, 
without fermions leaking out of the system, cf. Fig. [Tib). 

In an experimental realization, the hole fraction x 
is controlled by the boson fraction (also amenable to 
postselection). Moreover, tjj ~ TUf/2 and U nn /J ~ 
[2u~ x - t- 2 - l]u//2, with r = t b /tf and u v = U vv /U fb 
(y = /, 6), can be tuned by using Feshbach resonances 
and by modifying the lattices for bosons and fermions, 
either in depth or relative position. A candidate system 
is a mixture of Ytterbium isotopes [25M27I ] : The total an- 
gular momentum of Yb is just given by the nuclear spin 
J, not influencing the interparticle interaction. A mix- 
ture of two spin states of fermionic 173 Yb (1 = 5/2) with 
bosonic 168 Yb or 174 Yb (both I = 0) is described by three 
positive s-wave scattering lengths, with ajf/afb = 5.2 
(1.4), atb/afb = 6.5 (0.8) and a/6 = 2.0nm (7.3nm) for 
i68 Yd ( 174 Yb)_|20|- Also optical Feshbach resonances 
are available [27| ■ Using a far off-resonant optical po- 
tential allows one to create practically equal traps for 
bosonic and fermionic isotopes without fine-tuning. Con- 
sidering the Mott regime, according to band structure 
calculations, typical Hubbard interaction parameters for 
Yb will roughly be on the order of fO kHz. The effec- 
tive tunneling matrix element t will be about 200 times 
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smaller. Such low energy scales are challenging experi- 
mentally. However, they are not very specific to our pro- 
posal, but rather are a generic consequence of the fact 
that the t-J physics is inevitably the physics of superex- 
change processes. Therefore, the temperatures needed 
here are only moderately lower than those that would 
be needed in an experimental setup where — somehow — 
a Mott insulator of fermionic atoms can be doped with 
real holes in a controlled fashion. While in the latter 
case XJH /V 3> 1 would be required, in our proposal one 
has U ff /t f > 2{t/J) rs 6 (combine ru f ~ 2t/J and 
U ih ^> t b ), reducing J by a factor of 6 only. 

At hole doping x < 0(0. 01) and temperatures k-gT < 
J, the homogeneous t-J model (j^J is known to give rise 
to (quasi-) long-range antiferromagnetic Neel order. At 
larger doping x ~ 0(0. 1) and even lower temperatures, 
the model is conjectured to possess yet another ordered 
phase, namely the d x 2_ y 2-w&ve superconducting one ob- 
served also in the cuprates [H, [l(|. Superfluidity of c- 
fermions (that is in the t-J model) is connected to super- 
fluidity of both /-fermions and 6-bosons j28[; both have 
to be probed. The bosonic superfluidity is related to 
a (quasi) condensate at quasimomcntum k = 0, visible 
in time-of-flight absorption images. 2 . In order to verify 
e?-wave superfluidity of the fermions one can resort to 
different schemes for measuring quasiparticle excitations 
(as well as their interference) that have been proposed 
for cold atom systems for this purpose [TJ [29[ . 

In the cuprates, at optimal doping, superconducting 
behavior appears at temperatures that are at least one 
order of magnitude smaller than J [l(], [2l[ . It is challeng- 
ing to achieve these temperatures with cold atoms, and 
novel cooling techniques, as those proposed in Refs. [i~4j . 
have to be implemented. The system has to be divided 
into a low-entropy part (lacking low-lying excitations) 
and a part carrying most of the entropy (a metallic shell 
or a boson-gas reservoir). The latter is removed before 
the system is adiabatically transferred into the desired 
state. In our fermion-boson system, a suitable gapped 
low-entropy phase would be an insulator of two fermions 
and an integer number of bosons at each site, to be cre- 
ated by making the trap steep and the lattice deep. 

Finally, we would like to mention that neutralizing the 
trapping potential the way described here is also possible 
and interesting for bosons. In Refs. [16| it has been shown 

2 Bosons and "f and 4- fermions can be imaged separately by se- 
lective absorption or Stern-Gerlach separation. 



that the presence of a trap can fundamentally change 
the nature of the transition between an insulating and 
a superfluid quantum phase. The reason is not that the 
trapped system is too small to sharply distinguish be- 
tween the two phases. It is more subtle: Because of the 
inhomogeneity, the transition occurs locally by displac- 
ing the interface separating spatial regions of different 
bulk phases. However, this spatial interface might be a 
smooth crossover rather than a sharp transition. A direct 
consequence can be that the adiabatic passage between 
two quantum phases is greatly facilitated by the trap. 
A simple lattice model that would allow to study such 
phenomena is given by spinless hard-core bosons on a 
square lattice with nearest-neighbor repulsion U nn and 
hopping t. It features a transition from a checkerboard 
insulator at half filling and small t/U nn to a superfluid 
phase [3fJ . The model (being equivalent to an XXZ spin 
1/2 model) can be realized experimentally by creating a 
mixed Mott insulator having one boson of either species 
s = 1 or 2 on each site [l7j . Mapping b\ = b\ib i2 , the sys- 
tem is described by a single type of hard-core boson only, 
with U nn = (t\ + 4)/U 12 - Atl/U 2 2, t = 2t x t 2 /U u and 
residual trapping potential Wi — Vn — V^. Here U ss >, t s , 
and Vi s denote on-site interaction, tunneling, and traps 
for the species s = 1 and 2. Now Wi can be tuned con- 
stant without making particles leak out. Thus, a ramp 
of t/U nn at different rates can now be performed both 
with and without parabolic potential, and the degree of 
adiabaticity in passing the transition can be compared. 

We have proposed a robust implementation of a quan- 
tum simulator for the homogeneous t-J model with 
well controlled hole doping, using a sample of ultra- 
cold bosonic and fermionic atoms in an optical lattice. 
We believe that — once the necessary temperatures are 
realized — our scheme can serve to gain crucial insight 
into the physics of strongly correlated quantum matter. 
Moreover, realizing a homogeneous bosonic system allows 
one to investigate experimentally the role of inhomogene- 
ity when passing from one quantum phase to another. 
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